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Abstract 
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' The nonsinglet spin dependent structure function 172(2;, Q^) is studied at small x within the 

, double logarithmic approximation of perturbative QCD. Both positive and negative signature 

^\ ' contributions are considered, and we find a power-like growth in l/x. We discuss how our result 

, fits into the Wandzura-Wilczek relation. 

p 1' 

(D ! 1 Introduction 

X^^ , 

The evolution of the spin dependent structure function g\ is well known. In order to follow the 
^ ■ evolution at x ~ 0(1) one can use the original Altarelli-Parisi equation|l|] or the anomalous dimen- 

, sions calculated inp], ^ . In the region of very small x it is possible to use the double logarithmic 

approximation |5|| to examine the small-x behaviour oi g\. In contrast to this, the situation with 
other polarized structure functions, in particular is not so clear at the moment. 

The problem comes from the fact that has no parton interpretation. g2 describes the scat- 
tering on a transversely polarized nucleon, and in terms of the Wilson operator product expansion 
(OPE) it corresponds to the off-diagonal element of the density matrix (A^(A')|0|A^(A)) with the 
helicities A' 7^ A. 

Let us recapitulate a number of well known facts about the relation between the functions g\ 
and g2 (see, for example, the review [^): 

i) in the lowest approximation for the deep inelastic scattering (DIS) on a free (on mass shell) 
quark, the function 52(2;, Q^) = vanishes identically ( seeQ for a review); 

ii) in the limit m} <^ Q^, the "massless" quark may have only the longitudinal polarization; so 
gi{x,Q^)+g2{x,Q^)=0§; 
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iii) the Wandzura-Wilczek relation|8|] states, for the contribution of twist-2 operators, that 

1 

9i{x,Q^)+g2{x,Q^)= f-9i{y,Q') ; (l) 

J y 

X 

iv) for the spin dependent structure function 52 the next (twist-3) contribution is not suppressed, 
even for large Q. 

The presence of a twist-3 contribution ^2^^ {x, Q"^) in g2{x, Q^) can be seen in several ways. First 
of all, from an explicit one loop calculation ||To| , ^] it is known that g2 does not satisfy Eq.(l). 
Next, turning to the small-x region, it is again evident that the relation (1) cannot be valid for the 
whole functions gi and g2'- gi corresponds to an odd-signature amplitude, while g2 consists of even 
and odd signature pieces (see for example |12[). Therefore gi and g2 are expected to have quite 
different rr-behaviour near x = 0, i.e. they are proportional to different powers of at small x 
(see iSHl)- 

In earlier papers Q a simple evolution equation for the moments of twist-3 parts of 52 
(^2^^ (x, Q^)) function has been written, and the anomalous dimensions have been calculated for 



the scattering off a quark target. Subsequently, it has been realized |15| that the derivation in Ref. 



y, 14| was incorrect, because quark operators and quark-gluon operators of the same twist and 
quantum numbers mix with each other. In p!^ a complete basis of operators has been suggested, 
and the mixing matrix has been calculated in [|lq ]. Similar calculations have later been perform by 



several other authors [17, 19 1. Next in ref. the anomalous dimensions were obtained using 



another gauge and renormalization scheme. The authors confirm the result of ref. ||l^, but they 
disagree with [19|. 



The most significant departure from the standard leading- twist evolution equations ||2|, 14| is 



that the number of operators which contribute to the structure function g^\x, Q"^) is not fixed but 
increases with the moment index n. Nevertheless, in the two limits: 

i) number of colours Nc — > 00, or 

ii) the (x)-moment n ^ 00 

it has been shown that the quark-gluon operators decouple from the quark operator evolution 
equation [^], and the asymptotic behaviour of g2{x, Q"^) in the region 1 — x <C 1 was derived. 

In the region of small x <C 1 another new feature arises. It is known that at small x the 
strong ordering of transverse momenta is violated |22, 13, 23, |^. Instead of kf^ ^ ^ one has 



^ft ^ '^^^i-i,t (with Xi <C Xi-i). As a result, double log contributions of the form (a^ In 1/x)" 
appear which cannot be summed up in the framework of the conventional log evolution. For gi 
this point was discussed in detail in |Q, |5| . 

In this paper we study the small-x behaviour of the nonsinglet polarized structure function g2, 
in very much the same way as in |Q, ^ gi has been investigated. We will sum up all the leading 
pQCD double logarithms (DL) of the form a" In*^ l/xln™" with m < n and k + m = 2n — 1. 

The outline of the paper is as follows. First in Sect. 2 we will consider the non-singlet part of 
g2 in the first loop leading log approximation (LLA) and show that for the massive quark target 
(or for the //-meson in the QED case) in the small x region ^2^^ = 0. Section 3 contains a detailed 
dicsussion of the two loop DL contribution, and in section 4 we derive the DL evolution of g2 
at small x. In section 5 we analyse the signature content of g2, and in section 6 we discuss the 
connection of our results with the Wandzura-Wilczek relation (1). In the final section we present 
a brief summary. 
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Figure 1: The two Born diagrams 



2 Definitions and The One Loop Leading Logarithm 

We start with a few definitions and general remarks. The spin dependent antisymmetric part of 
the hadronic tensor ^W/^ of the DIS lepton-hadron amphtude has the form: 



~ (PQ) 



where M, Pf^, and denote the mass, the four momentum and the spin vector of the target, resp. 
As usual, (sp) = and = —1, and the spin-dependent structure functions are defined as the 
discontinuities of Ti and T2: 

91,2 = --^ImTi^2- (3) 
zvr 



Throughout our paper we work in the Feynman gauge. Note that with the help of the identity [12| 



(where the last two terms may be dropped due to the gauge invariance with respect to the photon) 
in (2) one can replace the vector in front of by —2xpaSfi. Then the spin structure corresponding 
to the function (72 takes the form {2xM) / {j)Q)e^^^°^l^ PaS(3g2- It is this form which appears naturally 
when we calculate the trace of the scattering amplitude. 

An important point that we have to discuss is the crossing symmetry of the amplitude. Begin- 
ning with the Born approximation we have to sum the two Feynman diagrams Fig. la and Fig. lb; 
they will be referred to as the s-channel and the u-channel graphs. The contribution of the u- 
channel graph may be obtained from the s-channel diagram just by permuting the photon indices 
H ^ V and by replacing Q — > —Q. The whole amplitude, given by the sum of Fig. la and Fig. lb 
contributions, should be symmetric under this combined ^ Q ^ —Q) transformation (this 
reflects the Bose statistics of two identical photons in t-channel). To provide such a property the 
function Ti should be odd, i.e. it must have negative signature. On the other hand, the crossing 
symmetry of the function T2, which is multiplied by the factor {2xM) / (pQ) ■ e^^^^^paSp is more 
complicated and, in fact, cannot simply be read off from (2). To illustrate the difficulty we simply 
note that, from the first part of the second term in (2) (the piece proportional to QaSp), we would 
conclude that T2 has the same signature as Ti. If, on the other hand, we make use of the identity 
(4) and neglect (by gauge invariance) the second term on the right hand side, the spin structure of 
T2 in (2) can be replaced by —2xpaSp. Now the same argument would lead to the conclusion that 
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Figure 2: The one-loop diagrams 



T2 has to be even, i.e. T2 has positive signature. In section 5 we will show that T2 has no definite 
signature but can be written as a sum of two pieces with opposite signature. 

To be definite let us consider the high energy scattering of a virtual photon on a heavy fermion 
with mass m. In order to obtain gi and §2 we are interested in the energy discontinuity. So in 
our calculations we will deal with the energy dicscontinuity, but in order to analyse the signature 
properties we will have to come back to the full amplitudes Ti and T2. 

From the tensor structure of (2) one sees that we have to consider a transversely polarized target. 
Namely, for a longitudinal polarization vector = the vector in front of T2 vanishes, and only 
Ti contributes. For a transverse polarization, on the other hand, both terms in (2) contribute, i.e. 
the result of such a calculation will be proportional to Ti + T2. We will therefore calculate the 
energy discontinuity of the scattering amplitude of the transverse polarized target which we will 
denote by g_\_. At the end we then subtract from the result the (already known Q) function gi. 

To begin with the Born approximation, we only mention that in the tree approximation the 
amplitude contributes only to the first term on the r.h.s. of (2). Therefore, the Born value of §2, 
92 Born, Vanishes. Contributions to the second term will appear only after the as corrections are 
taken into account. 



To study the one loop contributions it is convenient to use Sudakov |24| variables 



xp + aQ' + ktu] d'^k = ^dxdad^kt; Q 



2p'Q' = s ; Q^ = Q' - xp 



12 



V 



12 



; 



(5) 



and to separate the transverse and longitudinal parts of the gluon propagator. In the Feynman 
gauge the spin part of gluon propagator is: 



d 



PuQ'u + QuP'. 



{p'Q') 



(6) 



For the case of (72 it is better to deal with the transverse polarization of initial quark, as for the 
longitudinal quark only the function gi gives the leading contribution to the amplitude (2). 

So the density matrix of the target quark is: 



^{p + m){l - 75S-^), 

and the spin dependent part of the trace in Fig. 2a takes the form: 

1 



dapTr 



■jaip + m)j5S -ff3{k + m)'~iu{Q + k + m)jn{k + m) 



(7) 



(8) 
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Similarly, for the Fig. 2b we obtain 



ja{p + m)j5S^j^{Q + p + m)ji3{Q + k + m)^fj,{k + m) 



(9) 



where the quark momentum is qi = Q + p. 



It is easy to check that the transverse part of the gluon propagator d^is = g^p gives a vanishing 
leading log contribution to Eqs.(8) and (9). Indeed, due to the identity 



a/3 



(10) 



the expression (8) equals to zero. The situation with Eq.(9) is a bit more complicated. 

Note that as we deal with = the polarization vectors of the photons (//, u) in (4) have 
to be different. One has to be longitudinal, whereas the other one is transverse and orthogonal to 
s^. If we chose = 7^, j^, = 7]] we find, according to (10), that Tr^ = 0. If, on the other hand, 
7i/ = ^■nd 7(u = we arrive at: 



Tr" = -Tr 



{p - m)-f5S^'y^{Q +p - m){Q +p - m){Q + k + m)'y^{k + m) 



(11) 



At small X, in the quark propagator we simply put 1/A;^ 2± l/k^. 

In order to obtain the leading logarithm dk^/kf we have to keep in (11) only the longitudinal 
[k^ ~ xp) part of quark momentum and, due to the smallness of the quark mass ^ k^, one 
may keep just the first power of m. Then Eq.(ll) takes the form 



Tr^ ~ -(1 - x)mTr [p75sSi^Q(<3 + %J 



= 



(12) 



as in our LLA limit (5 + A;~(5 + a;p = (5' and pQQ' 
we neglect the contributions of the order of m?. 



xp^Q' + pQ''^ = —xni^Q' ~ 0; recall that 



Now let us consider the longitudinal part of the gluon propagator, i.e. the second term of (6). 
Keeping only the lowest power of mass the mass m, and taking into account that 



we obtain for the diagram shown in Fig. 2a: 



2m75S 



Tr"- = -mTr 



lbS^klv{Q + k)^^i.k 



= -mk'^Tr 



-m2{ks^)Tr ^^^AQ + k)^^k 



(13) 



Here one has two quark propagators (1/fe^)^ 1/^t • So in order to obtain the leading logarithm we 
have to keep in the numerator the terms proportional to A;|. After the averaging over the azimuthal 
angle of the vector kt the last term in (13) may be written as 



2m{ks^)Tr ^5lv{Q + k^)l^tk 



mk Tr 



+ mk Tr 



Now the first term cancels against the corresponding term in (13). In the last term, within the 
LLA, we may put k = xp, and finally obtain 



Tr" = +mxk^Tr 



(14) 
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For the case of Fig.2b the first longitudinal term daf3 = p'aQ'fs/ip'Q') gives zero (within the 
LLA). Thus we have 



2{p'Q' 



-Tr 



Q'{p + m)j5S-^{Q + P + m)p'{Q + k + m)'y^{k + m) 



(15) 



As before, in order to obtain the leading logarithm we put k = xp. Next it is useful to note that 
{Q+P + m)p'{Q + k + m) = 2{p'Q'){Q' + m). Therefore 



Tr" = -Tr 



Q'ip + m)75S-^7,,Q'7^(xp + m) 



-mTr 



Q'plbS IvQ'lii - {2pQ)mxTr 



mxTr 



(16) 



mTr 



X. 



The first term on the r.h.s. of (16) is equal to zero if 7^ = 7^. If on the other hand we take 7^ = jfi 
ilu = 7,}) we find 



-Tr 



Q'pi5S^luQ'l,^ = 2{pQ)Tr 75S^7.7)iQ' = -x{2pQ)Tr 



lul'ilP 



To obtain the last equality we have used the fact that the longitudinal photon polarization vector 

Jl _ 1 



■■{Q' + xp')^, i.e. 7/!||((5' + xp')fj_. Therefore in the last expression one can replace Q' by 

—xp. The factor 2pQ in Tr^ (15) cancels the rightmost quark propag ator l/{Q+pf ~ l/{2pQ). 

Finally, we have to add the contribution of the diagram 2c. After combining it with our results 
for Figs. 2a and b the whole one-loop LLA result takes the form 



im 

JpQ) 



■ {91 +92) 



MpQ) 



m 



d^Tr 



l^s^lf^luQ [(l + 0(x)]. 



(17) 



where Cq is the electric charge of the quark (7, Cp = (N^ — l)/2Nc, and is the QCD coupling 
constant. 

The well known result for the non-singlet function gi{x, Q^) in the one loop approximation and 
LLA small x limit is:: 



9i 



elCpas f dkf 



27r 



(18) 



Comparing Eqs.(17) and (18) we see that in this limit §2 = + 0{x), in agreement with the more 
precise computation [^] 

92{x,Q') 0. ^^xln , f , . (19) 



27r 



m'^x{l — x) 



3 Second loop 

Next we turn to the two loop contributions to §2- The diagrams are shown in Fig. 3. In this section 
we will show that only the first three diagrams in Figs. 3 a -c lead to a double logarithms, whereas 
in all other diagrams the double logarithm cancels. These cancellations are mainly due to the 
AFS-mechanism ^]: first one realizes that, in the approximation which leads to the double 
logarithms, it is enough to keep the nonzero value of the x variable only in those propagators which 
have the largest a. We then find for the integration of either xi or X2 that the propagator poles lie 
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Figure 4: Structure of the Sudakov vertex 



on the same side of the integration contour. Their sum adds up to zero, since closing the contour 
in the opposite half plane (which in this approximation has no propagator poles) gives zero. 

To see this in detail, let us go through the diagrams in Fig. 3. First we eliminate all self energy 
insertions to Fig. 2a or b: they never change the spin structure and, in the Feynman gauge, the self 
energy loop has only a single logarithmic divergence. Examples are shown in Figs.Sn and 3o. 

Next we consider the vertex correction in Fig. 3d. The absence of the double logarithm in the 
last step can be checked by a straightforward computation of the trace (analogous to (9), (15) for the 
case of Fig. 2b). Alternatively, one can use the " AFS" type arguments and consider the cancellation 
between the cuts passing through the quark lines with qi or g2- 

Next we consider the graphs in Figs.3e-h. In order to obtain the double log contribution we have 
to consider the domain where all the momenta are strongly ordered. Say, in Fig. 3g 1 ^ X2 3> xi, 
h2 ^ hi and I ail \a2\- Let us consider the loop with the largest values of the a component 
(the photon momentum fraction). This loop includes 3 quark propagators {qi,q2 and ki) and one 
gluon (og = ai — a2). Two propagators with the largest ~ ~ 1 are enough to provide the 
convergency of the integral over dxi. Thus we may neglect all other propagators (with the order of 
ai accuracy). For the l/qf and l/q2 propagators both poles lie in the same semiplane, and closing 
the xi contour in the opposite semiplane one gets zero result. In other words, the contributions of 
the l/ql and l/q2 poles cancel each other. This is analogous to the famous cancellation of the AFS 



pomeron cut contribution which was first discussed by Mandelstam |26]. This cancellation 
reflects the fact that the diagram shown in Fig.3g has a wrong space-time structure. In the Breit 
(photon rest) frame the formation time of the gluon with momentum g' = p — A;2, r' ~ l/{a'\/Q^), 
which is proportional to the inverse power of the a' ~ Q2 component r' ~ l/{a'\/Q^) is much larger 
than the formation time of the second gluon with q = k2 — ki, (r ~ l/(ai\/Q^), and also much 
larger than the s-channel quark time life, Tq ~ 1/ y/Q"^) (the separation between the points and 
V in Fig.l). On the other hand, along the upper quark line in Fig.3g we have quite the opposite, 
wrong time ordering > r > r'. By the same reason the graphs of Fig. 3 e,f graphs give no leading 
log contribution. In Fig.3e we have — ai = |afcj| ^ Oi'i- Then in the right hand side loop the two 
propagators with the largest a-components are §2 (og ~ 1) and /c2 = 92 — 9i ("2 ~ ai)- In both 
propagators the photon momentum fraction a goes in the same direction, i.e. we have both poles 
in the same X2-semiplane, and closing the contour of the X2 integration in the opposite semiplane 
one obtains zero. In the case of large a'^ ^> —ai (of course, we still have a'l <C 1 and — ai <C 1) 
we obtain the same zero result, due to the xi-integration in the left hand side loop. Finally, for the 
diagrams Fig.3f and h we consider the upper (A;i, gi, ^2, ^i) loop (without the quark k2): the two 
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propagators with the largest photon momentum fraction a are qi and q2, and we have the same 
cancehation as for the case of Fig.Sg. 

The diagrams in Figs.3i-k contain vertex corrections as subgraphs. The general structure of 
such a 'Sudakov' vertex is illustrated in Fig. 4. The double logarithmic (DL) contribution has 
been analysed first in [^J]. It comes from the region of ai ^> Og ^> 02, a;i <C <C X2 (for the 
case of ai 3> 02, xi <C X2). The corresponding kinematical domain for the graphs of Fig. 3 i,j,k is 
1 S> aq » a2i with q = k2 — ki in Fig.3i and q = qi—q2 in Fig.3j,k. Indeed, say, in the upper triangle 
of Fig.3i one may close the xi-contour around the 1/q'^ pole of the gluon propagator. Then the quark 
propagators 1/A;f ^ l/{x2aqs) (here xi = X2 — Xq ^ X2 as Xg <C X2) and l/qf ^ l/xgS. To obtain 
the leading logarithm one has to consider the longitudinal polarization of gluons d}j^p = Q'aP'is/iQ'p') 
(as we have seen in Sec. 2). In the trace corresponding to the quark loop it gives Q'k2Q' = X2S ■ Q' 
and p'q2p' = s ■ p' . In other words, we reproduce the original spin structure of the diagram Fig. 2b 
and obtain the factor X2 which cancels the l/x2 coming from the l/kf propagator. Thus the loop 
integral 

f dk'^Aandxa ^ / , x9\ f dai dxg , , 

/ " ' ' b [aqXqS + (A;i - k2ft) = / ^ ^ (20) 

takes the double logarithmic form. However in such a domain for the second loop we have the 
"AFS"-type cancellation again. In Fig.3i for the loop with 2 gluons and quarks k2 and q2 the two 
largest a (a^ and ^ 1) go in the same direction, leading to two poles in the same X2 semiplane. 
In the left loop of Fig.3j (which does not include the quark (72) the two largest a variables are Og^ 
and ag which flow in the same direction as well. For the case of Fig. 3k one has the same cancellation 
as in Fig.3g. In the lower triangle (with 3 quark and one gluon propagators) the two largest a are 
Qg^ ~ 1 and — 1. Here (in the vertex correction) the absence of a DL contribution may be 
treated as a well known cancellation of the leading logs between the virtual (Sudakov formfactor) 
and real "soft" gluon emission in the inclusive cross section; i.e. the cancellation between the cuts 
shown by the left- and right-dotted lines in Fig. 3 i,j,k. 

For the case of the vertices shown in Fig. 31, it is more convenient to fix the value of x' within 
the conventional DL region (1 ^ x' ^ xi) and then to observe the "AFS"-type cancellation in the 
upper loop {ki,qi,qg, but without the A:'-propagator) due to the two largest fraction of the target 
momentum Xg/ = x' — xi and Xg ^ l{qg = q2 — qi = XgP + OgQ' + qt) which go in the same direction. 

The vertex Fig. 3m does not give the double log at all. For x' <C 1 and a' <^ ai the virtuality 
{k' — ki)"^ of the quark propagator is controlled not by the value of aix' s but by k\^ ~ ais 0. Finally, 
in the case of the two photon vertex corrections (Fig.3q) one has to "cut" one of the vertices in 
order to have a small x. Then the central quark virtuality is ~ s, and again we loose the double 
logarithm in the second vertex, since the virtuality {qi — A;')^ ~ s of the next quark propagator is 
not controled by the gluon loop (a',x') variables. The same is true for the diagrams in Figs.3j and 
3r. 

Finally we return to the 'reducible' graphs of Fig.3a-c. They all have the property that they 
can be divided into two parts by cutting two t-channel quark lines. First of all we will check the 
cancellation of the double logarithms corresponding to the lower (below the crosses) part of the 
graphs in Fig. 3b and c, i.e. the cancellation at the level of the quark-quark scattering amplitude. 

■^If instead of the gluon qg one cuts the quark hues (p — k') and (fci — k') (as in Fig. 5a) we loose the logarithm in 
the upper loop (with gluon qg). Indeed, the quarks (p — k') and (fci — k') are now on mass shell, and thanks to the 
gauge invariance one may replace the polarization vector Q' (at the lower end of the gluon qg) by —qgijo-g — —kxtjcxx. 
This factor (fcit/ai) kills the normal Dlog in the loop. 
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Such a cancellation was first observed in QED |22| and was then discussed in |13] for the DL 
quark-quark amplitude in QCD. 

It is known that the DL contribution comes from the longitudinal spin part of the gluon 
k' propagator in the region of x ^ x' ^ 1 and |a'| ^ |ai|. The contour of the x' integration may 
be closed around the quark l/k^ pole, and then the contour of the next Xg integral (see Fig.5a,b) 
should be closed around the l/k^ pole (or one has to use the l/{kl + ie) propagator to take the 
imaginary part of the whole amplitude). This gives the "Sudakov"-type double logarithm 

similar to Eq.(20), but using the 5-function 6{kf) with ^3 ~ aix's + k'^^ for the d?k^ integral 
{k^t = k't — kit) rather than the sdx' integral. In terms of the s-channel intermediate states it 
means that we consider the cut of Fig. 3b, c diagrams shown in Fig.5a,b by the dash-dotted line; 
the upper part of Fig. 3b, c is shown in Fig. 5 by dashed lines. Note that in Fig. 5 we have crossed 
the upper quark line in order to plot both the graphs in the same form. The even signature 
(corresponding to 52 structure function) amplitude is invariant under the crossing transformation, 
and the only difference between Fig. 5a and 5b (or 3b and 3c) comes from the colour charge of the 
^3 quark line. The two t-channel quarks form a singlet colourless state and thus should have an 
opposite colour sign. Therefore, within the DL accuracy Fig. 3b and 3c, diagrams cancel each other. 

Note that such cancellation takes place not only in the case of Fig. 3b, c but for any DL quark- 
quark amplitude (in Fig. 5 the quark-quark amplitude is denoted by the dotted circle). This can- 
cellation, i.e. the absence of the non-ladder graphs in the DL even signature quark-quark QCD 



amplitude was proven and discussed in more detail in [13|. 



As a result of all these cancellations, the only DL two loop contribution to the even signature 
amplitude comes from the simple diagram in Fig. 3a with the two-quark t-channel intermediate 
state. However, the calculation of this pure ladder graph Fig. 3a is not so trivial, as we have to 
take care of the longitudinal parts of the quark momenta ki and /c2- For example, even in the DL 
kinematics the product of 4 vectors 2{kik2) = kf + k2 — Qg = kf + k2 7^ 2{kit ■ k2t)', here we put the 
gluon = on mass shell. That is why we will describe the algorithm of the computation in the 
next section. 
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Figure 6: Branches in the ladder graphs 



4 Evolution of Ladder Graphs 

As we have shown, in the two loop approximation only the diagrams Fig.3a-c give double logarithms. 
In this section we generalize our calculation of the ladder graph (Fig. 3a) to all orders (Fig. 6). In 
the next section we will then return to the non-ladder graphs in Fig.3b,c. 

We begin at the bottom part of Fig.6a and work our way upwards. After the first gluon rung 
the spin dependent part of the density matrix (7) takes the form 



a 



la 



p + m 



75 S la 



777.755 



(21) 



As before we keep only the first power of mass m <^ and therefore one may neglect the mass 

in all other propagators. Then we write the evolution kernel as 



with 



Cpots f dx dkf ^. 



^{ibs) = k^^sk = (fc^)75S - 2(sfc) • 75/c . 



(22) 



(23) 



The first term on the r.h.s. of (23) repeats the initial spin structure (75s) while the second one 
has the new matrix structure 75A:. When upper rungs are included, this second term generates 
a "new branch" (Fig.6b). Let us first follow this new branch. As there are only two external 
momenta (Q and k) available for this new term, the subsequent DL evolution with 77—1 rungs 
gives a contribution of the form 



2{s^k)-e^''''^k^Qp- f^'^-^Xk^Q) 



(24) 



where the tensor structure comes from the final trace, and the scalar function ^\k, Q) is equal 
(in (n-l)th order as) to the ladder part of the structure function gi{x' , Q^) (with x' = Q^/2{k ■ Q)). 
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Since the first rung in (22) and (23) has the same color factor CpOis/^TT as the evolution kernel 
in the ladder part of the longitudinal structure function (see below), we can absorb this rung into 
and obtain /^"^. Finally, from the integration over the azimuthal direction of the vector kt 
the product 2{s-^k)e^'^"^kaQp turns into kfe^^'^'^^ s-^Qj^. In this way we obtain the factor fcf which 
is needed for the logarithmic integral k^dk^/k'^ in the lowest cell of the ladder. As a result, after 
summation over n, the "new branch" resulting from the second term in (23) leads to a contribution 



of the form 



/(x,Q^ 




foH/2 



uj - fo{uj)/2 



CO 



ln:i + E(CF«.)"5?^(x,Q2) , (25) 



n=2 



where instead of the amplitude /q ^ in Q we now have put the even signature quark-quark ampli- 
tude: 



foiu) = Lo - \l uj"^ - as , (26) 

which describes the pure ladder DL evolution and does not contain the nonladder graphs. 

Let us now return to the "normal" branch in (23) (Fig. 7) with the structure 755-*-. At the 
next-to- lowest rung we simply repeat the argument in (23): we again find, in addition to the 
'old' structure, the beginning of a new one branch which, after multiplication with the two lower 
rungs and summation over n, again leads to a contribution (25). Moving up the ladder in Fig. 6, 
this procedure repeats itself at each rung, i.e. for each order as we obtain a contribution (25). 
Formally, this multiplicity of (25) is achieved by taking the derivative with respect to ag and then 
multiplying with a^: 

-c^s^g"- (27) 

oas 

where the minus sign reflects the sign of the second term in (23). 

Finally we have to complete the 'normal' branch which is obtained by following, at each rung, 
the tensor structure 755"*". As we have already said before, for each rung we collect the same 
coefficient as in the case of the longitudinal target polarization 0. Having reached the upper end 
of the ladder, we arrive at the "coefficient function" of the "normal branch" mj^s'^, which results 
from taking the trace at the upper end of the ladder. This trace 



Tr 



is by a factor 1/2 smaller than the trace Tr[j^s^'ji,Q'jfj] in the case of the function gi. The DLA 
result for the (ladder part of) spin dependent non-singlet structure function is given in (25). So 
the result for the "normal branch" is simply ^g^. 

*For the longitudinal case the analogue of (23) is k-y^pk = (fc^)75p — 2(pfc)75fc. In the leading-log kinematics we 
have 2(pA:) = as ~ kt , and in the last term we have to keep only the longitudinal part of fc ~ xp. Thus one may 
neglect this last {2{pk)^c^k ~ k^ ■ xj^p) term in comparison with the first {k^ ■ "fsp) one, and there is no 'new branch' 
in the longitudinal case. 
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Figure 7: Signature structure of the first loop 



Taking the sum of the ah contributions, we arrive at our final result for the ladder graphs in 
Fig.6: 



Since we have studied the transverse polarization of the target, our result represents the ladder- 
graph contribution to the sum gi + §2- Before we subtract gi, we first turn to the question of 
signature and the nonladder graphs. 

5 Signature and the Non-Ladder Diagrams 

So far we have generalized our two loop calculation only to the summation of ladder graphs. Now 
we have to consider the non-ladder graphs Figs. 3b and c and discuss their generalization to higher 
orders. In section 3 we have argued that the double logarithms of Figs. 3b and c cancel if we take 
the even signature combination of Figs. 3b, c and their s — u crossed counterparts, but add up if 
we consider the odd signature combination (Fig. 7). In this section we will analyse the signature 
content of the ladder graphs: we will find that a general ladder graph of Fig.6 contains both even 
and odd signature pieces. For the even signature contributions we then know that the nonladder 
contributions cancel, and only the odd signature pieces of the ladder graphs will receive extra 
contributions from the nonladder graphs. The summation of the nonplanar diagrams is done in the 
usual way of [|l3|,|, |. 

Let us therefore take a closer look at the signature properties of the ladder graphs. For this it 
is useful to discuss the amplitudes rather than the energy discontinuities. At the bottom of Fig.6 
we have the density matrix (21), mj^s. The contributions of Figs. 7a and b to the amplitudes Ti 
and T2 are proportional to: 




(28) 



1 



(29) 



2Qk + Q^ + k^ 



and 
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denominators we retain only the 2Qk term. The contributions (29) and (30) then take the form: 



Tr 

2Qk 



k\^5S - 2{sk)j5k)j,{Q + k)j^, (31) 



and 



2Qk 



Tr 



k\j5S-2{sk)j5khAQ-k)j^ (32) 



(in the last equation we have written the 7-matrices in the same order as in (31)). One immedi- 
ately sees that for the first piece (75s) two contributions with opposite signature emerge: the first 
one (proportional to r?'[75S7i/Q7^]) has odd signature, whereas the second one (proportional to 
Tr [75 S7jy/c7^]) has even signature. The second term (~ 2{sk)^^k) again leads to two contributions, 
but the second one vanishes (it has two identical k vectors, and the trace with the 75 matrix van- 
ishes). So we have only 2{sk)Tr[^^k'~^^Q'yf^\, and this structure has odd signature. Since this is the 
term which generates the 'new branch' we conclude that the 'new branch' is purely odd, whereas 
the 'old branch' contains both signatures. 

Next we have to relate this decomposition to the tensor structures appearing in (2). Let us 
first demonstrate explicitly that gi has only odd signature (cf. the discussion after (4)). For this 
we briefly return to the longitudinal polarization of the target. In footnote 3 we have argued that 
for Sn = Pn the second term in the decomposition (23), 2{pk)j5k, can be neglected in comparison 
with the first one, 75^. In the remaining term Trlj^pj^^lQ + k)j^] we neglect the k term, since 
in k = XkP + cxkQ' + kt the first piece gives zero contribution, and the second and the third one 
destroy the logarithmic divergence of the fct-integral. So, indeed, the longitudinal polarization is 
entirely odd signature, and gi has only odd signature. 

With this knowledge we return to the transverse target polarization and show how the varios 
pieces resulting from (31) contribute to the tensors in (2). The first term containing Q leads to 
the spin structure: Tr[j^s^^Q^^] = ^ie^vapQ'^s^ which contributes to both gi and g2- As said 
before, it has odd signature. The second (even signature) term has the spin structure of T2 (52): 
TrYi^s^i,k'yij\ = 'Ue^vaisk'^s^ , and in the DLA we approximate k^ XkP^ (cf. the discussion after 
(4)). Finally the spin structure of the last (odd signature) piece, 2{sk)Tr[y^k^vQl^i\- In order to 
have two logarithms we need a logarithmic integration k^dk^ / {k^)"^ . For the transverse polarization 

= we have already a kt factor from the prefactor (s/c), and in the trace we therefore have 

to retain the transverse component of k. Recall that this piece contributes to g± = gi + g2- On 
the other hand, for the longitudinal polarization = p^ we put (sk) = (pk) = afcs/2 = kf/2, and 
inside the trace we substitute k = x^p- This contribution is therefore smaller by the factor ~ xi^j 
than the 'normal' double logarithmic result for the longitudinal polarization. As a result, this last 
(odd signature) piece is negligable inside gi and contributes only to g2 ■ Combining all these results 
we arrive at the conclusion that gi has only odd signature contributions whereas 52 contains both 
even and odd signature pieces. 

To complete our discussion we have to consider diagrams with one more rung (Fig.8). Prom the 
density matrix at the bottom we now have not only the 'old' matrix 75s but also a term of the form 
(sZ)75l. For the former we can repeat all previous steps. For the latter, however, it is nontrivial 
to show that, at the double log level, the term {sl)Tr[k^5lk^,^k^fj] does not contribute to the even 
signature case (otherwise not only g2, but also gi would receive an even signature contribution). 
The even signature trace term is written in the form 

{sl)Tr[k-/Sluhn] = k^{sl)Tr[-idlvkln]- (33) 
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Figure 8: Signature structure of the second loop 



One possibility for obtaining the integral dlf/lf is to keep the transverse components in both 
Z-f actors. This gives: 



(34) 



The other possibility is to retain, in the trace, the large (x; S> x^) longitudinal component of l^. 
Instead of the rhs of (34) we obtain: 



k^{sl)Tr[^^xip^„ka^i\. 



(35) 



At first sight this expression does not seem to have a double logarithm from the It integral. On 
the other hand, it is larger than (34) by the ratio xi/xk, and we have to remember that there are 



corrections from the propagators 1/k . Putting k = a^x^s + kf with 

1 



xis 



we expand: 



I Xk ik-l)t 



1 + 



2(/ct/t) Xk 



XI 



+ 



(36) 



Combining the correction term on the rhs with (35) we obtain the double logarithmic contribution: 



k\sk 



2{ktlt 



■Tr[-ir>xfrivk-if, 



Tr[y5XkP'yuS-ftj] 



(37) 



which exactly cancels the contribution from (34). This shows that, in fact, also at the two loop 
level the structure (5^)75/ leads to odd signature only. 

From this discussion the following general pattern emerges. As we have shown before, in a 
general ladder diagram of Fig.6 each rung reproduces the tensor structure 75s and creates, in 
addition, the beginning a new branch. Follwing the first branch, we have a mixture of even and 
odd signature, whereas the new branch belongs to odd signature only. As a general result, during 
the evolution along any of these branches signature is conserved. For the first branch, the final trace 
at the top of the ladder separates even and odd signature parts. In order to include the nonladder 
graphs in the correct way we have to discuss even and odd signature branches separately. 

Let us first follow the odd signature part of the branch (755-*-): at each step it creates a new 
branch of pure odd signature with the spin structure 75^. For any odd signature branch in Fig.7 
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(say, with the structure 755-*- up to the j-th rung, and with the structure 75^; above the j-th rung) we 
have the well-known infrared evolution equation which takes into account both the ladder graphs 
and the double logarithmic non-ladder bremsstrahlung contributions. For each such branch we 
have the same result as for gi: 



Here \lo) denotes the quark scattering amplitude: 



f(-\u) = u - - + /^^. (39) 



The non-ladder diagrams are contained in the function 

■ d_ 
du 

where 



= 47riVaf A In {e^"'^D,{z)) , (40) 



^ = "/V^ (41) 

(for further details sec the Appendix). Once we reach the top of the ladder, we 'close' with the 
trace. The odd signature part is induced by the Q part. 

The correct counting of these odd signature branches is done in the following way. New branches 
are created only at ladder rungs. Non-ladder gluons never change the spin structure. We therefore 
have to count the number of rungs, and we do this by formally distinguishing between (the 
strong coupling multiplying a gluon rung) and (the strong coupling multiplying a nonladder 
gluon). This procedure turns the in the second term in the square root expression (39) into a^, 
whereas the ag in the third term of this square root, as well as in (40) and (41) becomes . As 
in (27), the number of new odd signature branches is obtained by applying the derivative operator 
— onto gi. In addition we have a contribution gi from the main 75s branch. So altogether 
we find: 

5l''=(l-«f^)5l|a.=«f=a.- (42) 

This generalizes the odd signature part of (28), in that in includes the nonladder diagrams. 

Next the even signature branch: it is contained only in the tensor structure 755^, and it is 
projected out by the k part of the final trace. Its coefficient is by a factor 2 smaller than in the 
usual gi calculation (since for the longitudinal photon one can replace Q by —2xp). For this branch 
nonladder graphs cancel, so our result equals: 



(43) 



(cf.(25)). 



Finally, in the one loop contribution we have to add also the nonladder graph. In section 3 we 
have argued that in the diagrams of Fig.3d the two cuts across qi and q2 cancel against each other. 
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However, in the first loop (Fig.2b,c), where the cut of the quark qi line corresponds to x = 1 and 
where we have the term fyy^s in the initial quark density matrix, the usual cancellation between 
the real and virtual DL contributions does not work. Hence within the DLA we have to take care 
about the non-ladder diagrams Fig. 2b, c. They give a contribution of the form: 

(44) 

Now we can add all contributions, (42), (43), and (44). Since this result for the transverse 
polarization is proportional to the sum g\ +52; we have to subtract the known result for g\. We 
obtain: 

(72(x, Q2) = -^/(x, Q2) _ (^LA_g^ Li=a^=«. + (45) 

Eq.(45) is our final result for g2. It differs from the one obtained in |27| in several respects: in [27| 
crossing properties and non-ladder graphs had not been analysed. As a result, the even signature 
contribution was omitted, and the second term did not correctly contain the nonladder graphs. 
Finally, the last term representing the one loop contribution was missing. 



6 The Wandzura-Wilczek Relation 

Let us finally comment on the Wandzura-Wilczek (WW) relation (1) and see how our results 
fit into this equation. Since the WW-relation was proven to be valid only for the twist-2 part 
of 52 but not for the whole structure function g^^ we do not expect that our DLA result for 
the small-x behavior satisfies this identity. It is instructive to illustrate the role of the twist-3 
part in lowest order perturbation theory. For small ^ 1 both functions g\ and g^ behave as 
const+0(as In (5^) + 0(as In i), but for 52 the "const" vanishes, and the term 0(as) goes to zero 
proportional to x. Therefore, at the order a^, the integral on the r.h.s. of eq.(l) contains one 
extra Inl/x, which on the Ihs can be contained only in the twist-2 part of 52- But since the full 
92 (6q-(19)) has no such a logarithm at this order of the coupling a^, this logarithm in the twist-2 
part must be cancelled by a similar logarithm in the twist-3 part of 52- This illustrates that in the 
small-x region the twist-2 and twist-3 contributions of 52 cannot be separated in an easy way. 

The systematic construction of the twist-2 and twist-3 operators has, for example, been de- 



scribed in |16|. The n- moment of the spin dependent structure function is given by the expectation 
value of the operator: 

^=(02) Q'^,,-Q'^A^l^l.'D^,...D^^i,\ , (46) 

where the matrix 70- is "directed" in accordance with the quark spin vector s^. In order to separate 
the twist-2 and twist-3 components, we have to symmetrize. The most symmetric operator 

^1 = Q',.-Q',r. Skm.,...,m.} [V^TsT.^m.-^m^] (47) 

gives the leading, twist-2 contribution, while the 

/ 2i \ " 
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corresponds to twist-3 (in eqs. (47) and (48) the symbols Sj^^^,,,^^^} and ■^[a,tj.i] denote symmetriza- 
tion and antisymmetrization with respect to the indices and [a, /xi], respectively.) In 

the case of longitudinal polarization of the target quark the vector has the same longitudinal 
direction as the derivative i^^. Therefore the operator i?2 does not contribute to the longitudi- 
nal polarization, and twist-3 does not appear. This situation changes when one deals with the 
transverse spin vector s^. Now the twist-2 contribution is given by the operator Ri, e.g. for n = 1: 



[h\Ri\h) r.Q'^(h\i^^^-f^{xp^)^h) + Q'^ih\i^-f^^^k^Mh) , (49) 
and i?2 becomes nonzero: 

{h\R2\h) ^Q'^(h\i;^5laixp^)i^\h) - Q'^{h\i;j5j^k^^l;\h) . (50) 
Now it is easy to see that the Wandzura-Wilczek relation is valid only for the twist-2 part of the 



structure function. It follows just from the relativistic invariance [16|; both sides are the matrix 
elements of the different components of the same twist-2 tensor operator. 

The results of our calculation show that at small x the WW-relation looks rather artificial. In 
particular, the twist-3 contribution is by no means small. For the first-loop we have shown that 
the twist-3 contribution adds an extra term to g2, which is even more singular at x ^ 0. For the 
case of the one loop 0{as) contribution the situation (and the behaviour of the moments of twist-2 
and twist-3 operators) has been discussed in detail in In the context of the MIT-bag model, a 
large twist-3 contribution to 52(2;) has been emphasized in |28]. As a result, we conclude that the 
WW relation cannot be used to estimate the full function §2 at small x. 



7 Conclusion 

In this paper we have calculated the small x asymptotics of the spin dependent structure function 
g2 in the double log approximation of perturbative QCD. We started with the one and two loop 
approximations, and we then succeeded to perform the sum over all orders in a^. The asymptotic 
small X behaviour of §2 is governed by the rightmost singularity in the w-plane at w = liJo — V^Cpo^ 
(~ 0.4 at as = 0.2). This is the square root singularity of gg-amplitude fo{uj). As a result, at x — > 
the structure function g2 behaves as 

52 « i^j (51) 

with the anomalous dimension 7 = a;o/2 . 

There is an interesting difference between our transverse distributions g±{x,Q'^) (or g2{x,Q'^)) 
and the structure function /ii(x,Q^) which describes the evolution of the operator (j^Q's^)- The 
small x behaviour of hi{x, Q"^) was studied in |29|. In the case of g2{x, Q^) the double logarithmic 



ladder evolution has been found to take place near the singularity at j = 0, and g2{x,Q ) ~ x~ , 
with A = Wo ~ 0{y/a^) > 0. In hi, on the other hand, the situation is quite different. Namely 
due to the fact that after the interaction with the s-channel gluon the spin structure vanishes 
{daffJalbQ' S'^7f3 = —^{Qs'^)l5 = 0), the leading ~ (l/x)*^ part of the function hi has no DL 
evolution at all, i.e. the singularity at j = in the complex angular momentum plane does not 
contribute. Only the singularity at j = —1 leads to the DL contribution, and the small-x behavior 
of hi is suppressed by one power of x compared to g2- /ii(x,Q^) = x/(asln^ 1/x, as In 1/x In Q^). 
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Let us recall that in this paper, as a first step, we have calculated g2 only for deep inelastic scat- 
tering on a quark target, i.e. only the evolution of the 52 structure function within the perturbative 
QCD. Nothing has been said about the initial distribution of (72 coming from the confinement re- 
gion. On the other hand, in ref.[^] it has been argued that, within the parton model, the intrinsic 
transverse momenta of the quarks (inside a proton) should lead to a non-zero value of 52 which 
satisfies the Wandzura-Wilczek relation (1). However, in order to reproduce the relation (1) one 
has to use very large transverse momenta {k^ ~ Q^; as the integral over corresponds in ||30|] to 
the integration over dy in eq.(l) with dy = dk^ /2mu). Within the parton model such large values 
of intrinsic momenta kt look somewhat unnatural. 



8 Appendix 

Here we will review the structure of the odd signature quark-quark amplitude and show where the 
QCD coupling (a^ = /At:) corresponds to the ladder, and where it belongs to the non-ladder 
contributions. We have to separate these contribution in order to take the derivative with respect 
to the ladder coupling in the final expresion (45). 

The well known equation for the odd signature quark-quark amplitude in the uo representation 
has the form[^ 

/o M-^^- + ^(/o Hj -^n^jih M (1) 

where the first two terms in the r.h.s. generate the ladder contribution while the last one corresponds 
to the non-ladder graphs. 

This last term on the rhs is due to the signature changing contributions of gluon bremsstrahlung 
which lead us to define an even-signature quark quark scattering amplitude f^\ijj) with color octet 
quantum number in the t-channel. This amplitude has the infrared evolution equation 

Its solution follows from the discussion given in [13|: using the transformation 

i'-\u^)=Nglj^^\nu{u) (3) 
the Riccati equation (3) is equivalent to the linear differential equation 



d'^u du 1 qi 

— K-z -4i^n = (4) 

dz^ dz 2iV2 gj^^ ^ ' 



where 



After a trivial transformation this equation is solved by a parabolic cylinder function. As a 
result, /g^"* has the form: 

i^\u:) = Ng%,^ln {e^'/'D,{z)) = Ng%,^ In (|^°^ dtr'-^e-^'-^y^^ (6) 
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with 

^ 9l 



2A^' 9nl 

With this we return to (1) and obtain for fQ~': 



(-). 



(the minus sign in front of the square root follows from the requirement that, for large co, the 
solution has to match the Born approximation). 



Acknowledgement: One of us (M.G.R.) gratefully acknowledges the hospitality of the DESY 
Theory Group. 

References 

[1] G.Altarclli and G.Parisi, Nucl.Phys. B126 (1977) 298. 
[2] M.A.Ahmed and E.G.Ross, Nucl.Phys. Bill (1976) 441. 
[3] R.Mertig and W.L.Van Neerven, Z. Phys. C70 (1996) 637. 

[4] J.Bartels,B.I.Ermolaev and M.G.Ryskin, DESY-95-124, Z. Phys. C70 (1996) 273. 

[5] J.Bartels, B.I.Ermolacv and M.G.Ryskin, DESY-96-025, Z.Phys. C7 

[6] M.Anselmino, A.V.Efremov and E.Leader, Phys.Rep. 261 (1995) 1. 

[7] B.L.Ioffe V.A.Khoze and L.N.Lipatov, "Hard Processes", North-Holland, 
Amsterdam, 1984. 

[8] S.Wandzura and F.Wilczek, Phys.Lett. B72 (1977)) 195. 

[9] A.J.G.Hey and J.E. Mandula, Phys.Rev. D5 (1972) 2610. 
[10] G.Altarclli, B.Lampc, P.Nason, G.Ridolfi, Phys.Lett. B334 (1994) 187. 
[11] A.Harindranath and Wei-Min Zhang, Phys.Lett. B408 (1997) 347. 
[12] R.L.Heimann, Nucl.Phys. B64 (1973) 429. 
[13] R.Kirschner and L.N.Lipatov, Nucl.Phys. B213 (1983) 122. 
[14] J.Kodaira et al., Phys. Rev. D20 (1979) 627; Nucl. Phys. B159 (1979) 99. 
[15] E.V.Shuryak and A.LVainshtein, Nucl. Phys. B199 (1982) 951; Nucl. Phys. B201 (1982) 141 



20 



[16] A.P.Bukhvostov, E.A.Kuraev and L.N.Lipatov, Sov.Phys.JETP 60 (1984) 22. 

[17] P.G.Ratcliffe, Nucl. Phys. B264 (1986) 493. 

[18] I.I.Balitsky and V.M.Braun, Nucl. Phys. B311 (1989) 541. 

[19] X.Ji and C.Chou, Phys. Rev. D42 (1990) 3637. 

[20] J.Kodaira, Y.Yasui and T.Uematsu,Phys.Lett. B344 (1995) 348; J.Kodaira, Y.Yasui, K.Tanaka 
and T.Uematsu, Phys.Lett. B387 (1996) 855. 

[21] A.Ah, V.M.Braun and G.Hiller, Phys. Lett. B266 (1991) 117. 

[22] V.N.Gribov, V.G.Gorshkov, G.V.Frolov and L.N.Lipatov, Sov.J.Nucl.Phys. 6 (1967) 95. 

[23] B.I.Ermolaev, S.I.Manayenkov and M.G.Ryskin, Z.Phys. C69 (1996) 259. 

[24] V.V. Sudakov, Sov.Phys.JETP 3 (1956) 65. 

[25] D.Amati, S.Fubini, A.Stangenini, Nuovo Cim. 26 (1962) 826. 

[26] S.Mandelstam, Nuovo Cim. 30 (1963) 1113, 1127. 

[27] B.I.Ermolaev and S.I.Troyan, |hep-ph/9703384 . 

[28] R.L.Jaffe and X.Ji, Phys. Rev. D43 (1991) 724. 

[29] R.Kirschner et al., Z.Phys. C74 (1997) 501. 

[30] J.D.Jackson, G.G.Ross and R.G.Roberts, Phys.Lett. B226 (1989) 159; 
R.G.Roberts and G.G.Ross, Phys. Lett. B373 (1996) 235. 



21 



